Oblique launching of optical surface waves by a subwavelength slit 
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The electromagnetic field on the metal surface launched by a subwavelength slit is analytically 
studied, for the case when the fundamental mode inside the slit has a wavevector component along 
the slit axis (conical mount). Both near- field and far- field regions are discussed, and the role of 
surface plasmon-polaritons and Norton waves is revealed. It is shown that the distance from the slit 
at which Norton waves are more intense than surface plasmons decrease with parallel wavevector. 
Additionally, it is found that the s-polarization component, while present for any non-zero parallel 
wavevector, only weakly contributes to the NWs. 
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I. INTRODUCTION 

Launching surface plasmon-polaritons (SPP) along 
metal surfaces has recently attracted a lot of inter- 
est for its possible application in integrated optical 
devices. One of the most common configurations 
used employs systems with translational symmetry 
in one direction, as a subwavelength slilP^ or a line 
defect^ ^. 

Very recently, several studies have been devoted to 
finding simple analytical models for the electromag- 
netic (EM) field radiated by a single slit, which pro- 
vides insight into the relevant physical pro cesses. ^^^^^ 
It has been found that the field at the surface presents 
a rich behavior as a function of both frequency and 
distance to the slit. Typically, at distances less that 
2-3 wavelengths, the field presents a complex spatial 
dependence^^ (which is sometimes phenomenologically 
described as composed of a SPP plus a "creeping" or 
"quasi-cilyndrical" wave, CW^. Atlar ger distances, 
there is an intermediate regime where the EM field 
is dominated by the SPP contribution, and a long- 
distance regime where the SPP has become negligi- 
ble and the field is that of a Norton wave (NW). All 
the previously cited works focused on the case of field 
propagation perpendicular to the slit, and did not ad- 
dressed the situation of non-normal incidence of light 
into the defect. This is a serious deficiency as SPP 
scattering effects are easier to detect when the SPP 
is launched obliquely, i.e., when its wavevector has a 
component along the slit axis.l^HI^ 

In this paper we fill this gap, presenting an anal- 
ysis of the EM fields on the metal surface obliquely 
launched by a subwavelength one-dimensional aper- 




FIG. 1: (color online). The geometry of the studied sys- 
tem. A plane electromagnetic wave impinges onto a sub- 
wavelength slit, placed in an optically thick film. The angle 
of incidence is arbitrary provided the electric field points 
perpendicular to the slit axis. The slit aperture at the 
transmission region generates a field at the surfaces, prop- 
agating non-perpendicular ly to the slit axis. The solid cir- 
cle represents a slice of the light cone for a fixed frequency, 
while the dashed one is the slice of the "SPP cone". Ar- 
rows inside the slit indicate the propagation direction of 
the fundamental mode. 

ture. Such a launching could be realized by illumi- 
nating a subwavelength slit in a thick metal film by a 
plane wave in the conical mount, where the incident 
wavevector and the electric field have a non-zero com- 
ponents along the slit axis and across the slit, respec- 
tively (see Fig. [T]). From now on, we will denote this 
configuration as "oblique incidence". The subwave- 
length slit "filters" the EM field inside it, so that only 
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a capacitor-like fundamental mode transfers the energy 
to the outgoing face of the film, all other modes inside 
the slit being exponentially suppressed. For oblique 
incidence, the /c- vector of the fundamental mode has 
a finite projection along the slit axis and therefore the 
field emerging from the outgoing aperture gains both 
p- and s-polarization components. Here we do not con- 
sider the part of the problem related to the transmis- 
sion and refiection efficiencies, but concentrate on the 
electric field pattern on the outgoing face of the film. 
We thus only need to know the Green's function of our 
system and the field at the outgoing face of the slit. For 
this we apply the mathematical methods described in 
Ref. 19, which have been previously used for obtaining 
the field at t he metal surface radiated both by a single 
slit for ky = (pEEel ^nd by a subwavelength hol^. 

Anticipating things, we would like to stress that the 
situation where ky ^ 0, presents two main differences 
with the case ky = 0: {i) the transmitted electric 
fields have components parallel to the slit axis and (n) 
the propagation length of the excited SPP along the 
direction across the slit diminishes approximately as 
Lspp oc cos Oi. The latter fact favors bringing closer 
to the slit the region where Norton waves dominate. 

II. ANALYTICAL SOLUTION 

Let us consider a plane monochromatic wave inci- 
dent onto a thick metallic film with a subwavelength 
slit centered at X = 0. The wave vector of the incident 
wave is k^, forming an angle Oi with respect to the 
Oy axis, and the wave is polarized perpendicularly to 
the slit (see Fig. [T]). If the incoming dielectric media 
has the dielectric permittivity £i, ky = -^/iik^ sin Oi, 
where k^ = 27r/A. Due to the translational symme- 
try of the system, the ^-component of the wavevector 
is conserved in the transmission process. This places 
constrains onto the angular interval where transmitted 
radiation can be found. If the outgoing half-space has 
dielectric permittivity £2, the range of the allowable 
angles is given by ky < k^^^Je^-, or sin^^ < ^Je2|s\ 
(otherwise, the field has evanescent character in the 
outgoing region). For example, if the field is incident 
from the glass substrate with e\ = 2.25 and the out- 
going medium is vacuum, the permissible angles are 
Oi < 41.81°. 

A. General analytical expression for the field 

Let us start with a note on the notation used: 
throughout the paper, distances in lower case let- 
ters are expressed in dimensionless units, defined as 
X = k^X, z = k^Z, and the dimensionless wavevector 
components are denoted as Qx^y = kx^y/k^ (so that. 



for instance, the light cone in vacuum corresponds to 

9 = 1). 

Starting from Lippmann-Schwinger integral equa- 
tion and taking into account that the field in 
the slit points along the x-direction E(x, y, z) 
exEx{x, z)e^^yy , the expression for the transmitted 
field simplifies to (see Appendix for details) 

E(x,z) = C j dx'G{x-x',z)Ex{x',z = -6), (1) 

where C = iyjsm — s/k^. The dielectric constants of 
dielectric in the outgoing region and the metal are e = 
£2 and respectively, and 5 is the skin depth in the 
metal. The integration in x' is performed across the 
slit area. The cyclic dependency upon the coordinate 
y, (X e^^yy is omitted here and in what follows. G(x, z) 
is the x-column of the Green's Dyadic, G = Gbx (with 
ex being unitary vector along Ox axis), whose angular 
representation for the case of arbitrary ky reads 

G^(x,z) = ^ / ^^t^ I IwA e^(^^-+^^^), 

G%x,z) = ^ / 4^t^ l-llqy] e^^^^-+^^^), 

(2) 

where G = G^ + G^ indices "p" and "5" stay for 
corresponding polarizations and tp and ts are Fresnel 
transmission coefficients for the metal-dielectric inter- 
face, given by 

t'= ^^"^ , t^ = ./^ ^^"^^ , (3) 

with q'^ = ql ^ q^, q^ = s - q^ , q^m = \/ Sm - q^ - 
The branches of qz, qzm must be chosen in accor- 
dance with the radiation conditions lni{qz,qzm) > 0. 
It should be noted that for qy = ^ the integrals trans- 
form to the case of in-plane launching (see Refs jl5|16|) 
and s-components of the fields vanish. 

For narrow (subwavelength) slits, the field inside 
the slit can be taken as independent of x and, there- 
fore, the outgoing slit aperture is equivalent to the ef- 
fective two-dimensional electric dipole located on the 
metal surfacc"*^^: E(x, z) = G{x, z)p^^^ , where p^^^ = 
aCEx{0,z = —5) and a is the width of the slit (in 
dimensionless units). 

The numerical computation of the integrals in 
Eq. ([2| are notoriously difficult, due to the simultane- 
ous presence of poles, branch cuts, and strongly oscil- 
latory factors. Using a special mathematical treatment 
based on the steepest descent method, an accurate 
analytical representation of the asymptotic behavior of 
the field is possible. We present here the final result for 
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the field at the metal surface z = 0; the mathematical 
details can be found in the Appendix. Introducing the 
following notation 



G{x,z = 0) = ^g(^), 



(4) 



we have in the region qox ^ 1: 

g(x) « iTrCpe^'^f^erfcf-is^Vft^) + e*«"^,/^^- 

\ qoX Sp 



2C 



^ +2\/2e-'T(f- +fP) 



(5) 



In this equation qo = yje — Qy presents the inverse 
spatial period along the x-axis of the algebraically- 
decaying terms, = Y^^^^rrzTl^^^^m) is the mod- 
ulus of the in-plane component of SPP momentum 

The fac- 



Qy IS its x-component. 



and Qxp 

t Of Sp = e~'^^^^y^qxp/qo — 1, appearing in the argu- 
ment of the complementary error function, erfc, is the 
position of the pole in the complex plane where the 
steepest-descent integration is made (see Appendix). 
It has an important significance, being responsible for 
the asymptotic expansion of the error function, whose 
argument is the the square root of the so called numeri- 
cal distance introduced by Sommerfeld. Another prop- 
erty of Sp is that quantify the distance in complex 
(7— space between the SPP pole and the branch-point 
placed at g = ^/e (i.e. at Qz =0). 

The terms in Eq. ([5| containing C„ come from the 
singular part of the integrals in Eq. pi) appearing for 
p-polarization. Cp are the residues at Qx = Qpx- 
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^ I QxpQzp \ 



(6) 



X ^ 1/ y^r^-^), similarly to what occurred in the case 

of Qy = G^^. Notice, nevertheless, that as Qy increases 
this formal condition is fulfilled for larger values of x, 
which explains that, at a fixed distance, the relative 
error in the field increases with Qy. To characterize 
the relative error, we have introduced the following 
function Ag^ = \{ga -9T'^)l9T'^\ with a = x,y,z, 
where g^'^^ correspond to precise numeric calculations 
and ga is given by Eqs. ([5|-([7|). Except for qy ^ 1 
this error is not very sensitive to the value of e^, and, 
therefore, almost independent of wavelength (from the 
optical region to longer wavelengths). We have checked 
that the relative error does not exceed a few percents 
for distances as small as X ~ A/10, and already for X 
of order of a wavelength the error has reduced down 
to 0.1% Sit qy=0 (and - 0.5% at Qy ^ 1), see Fig.|2j 
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FIG. 2: (color online). The relative error Aqx as a func- 
tion of Qy in the case of a gold surface. Continuous curves 
correspond to A = 700nm, while discontinuous ones are for 
A = 540nm. 



where q^p = s/^e ^Sm- 

Finally, f^'^ are contributions to the algebraically- 
decaying term, that in this ID geometry go like ^ 
1/x^/^, which dominate the far-field region at the sur- 
face (unless the system is completely absorptionless) 



(7) 



It is remarkable that the region for validity of the so- 
lution Q-Q is much less restrictive than q^x ^ 1 (i.e. 




B. Perfect electric conductor limit 

Before analyzing the case of a real conductor, let 
us consider the limiting case of a perfect electric con- 
ductor (PEC), characterized by Sm —00. Then the 
electric field at the surface is perpendicular to it and, 
for a very thin slit, can be analytically calculated using 
directly Eqs. ([lj,([2|: 



E 



PEC 



(x,0) = -e,^aHi'\qo\x\)E^^^{0,0), (8) 



with h[^^ being Hankel function of first order. 

The PEC limit of the asymptotic expansion has to 
be taken with care, as the Green's Dyadic tends to 
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zero as G ~ This is natural, since the dipole 

placed on the metal interface and oriented along it, 
cannot radiate due to cancellation of the field by the 
image dipole. However, the effective dipole of the slit 
diverges as p^^^ ^ V^^, so that the product G p^^^ re- 
mains finite. If we substitute Eqs. Q-Q into Eq. ([T]) 
and perform the PEC limit, we arrive at the asymp- 
totic expansion for Eq. ([8|. This expansion consists of 
just one term, where the Hankel function is replaced by 
its asymptotic term H[^\qo\x\) = —^J2il(i\q<^x)e^'^^^. 
Thus, our asymptotic expansion recovers the PEC re- 
sult, up to terms of order 0(x~^/^). 



C. Far- field asymptotic 

In case of long distances, or, more precisely, when 
xgol^pl 1, we can obtain a simplified expression from 
Eq. (|5| (exact up to 0(x-^/2)): 



g(x) = g5Pp(^) + gw(^), 



(9) 



where gs-pp is the contribution from the SPP pole 
(arising from the first term in the asymptotic expan- 
sion of the complementary error function) 



g5Pp(^) 



(10) 



and gAriy(x) = g^^(x) + g|^^(x) is an algebraically- 
decaying term, with contribution from both polariza- 
tions 



gArH^(^) 



/7re 



a =p,s. 



(11) 



This field component presents the two-dimensional 
analog of the Norton waves (NW) discovered theoret- 
ically almost a century ago by Norton^^, when ana- 
lyzing the radiation of radio waves by a point dipole 
placed over the Earth surface (represented by a lossy 
dielectric). 

We would like to stress that the approximation (|9|- 
( pTj ) is not applicable for PECs since, in this case, 
Sp = and the condition for validity of the asymp- 
totic expansion is not fulfilled for any distance. 

In the remaining part of the article we analyze the 
dependency of the fields upon qy. We concentrate on 
the case of a vacuum-gold interface setting £ = 1, 
= ^Au- The effects generated by changing e were 
considered in Ref. |16| 



III. ILLUSTRATIVE RESULTS AND THEIR 
DISCUSSION 

Before presenting the dependence of the field pat- 
tern on qy^ for completeness and in order to make the 



comparative analysis easier, we briefiy review the case 
qy = (extensively studied in [T4HT7|) . Figure [s] shows 
a representative case (A = 540 nm and qy = 0) for the 
dependence on distance to the slit of the electric field 
at the metal surface. The fields for both vacuum-gold 
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FIG. 3: (color online). The dependency upon distance from 
the slit of the electric field modulus at the gold surface. 
The wavelength is A = 540 nm and Qy = 0. Together 
with the total field, both SPP and NW contributions, and 
the PEC case are also plotted. The insets show the real 
part of the z-component of the total, SPP and NW fields, 
in different spatial regions. All fields are normalized to 
E{X = 0.1X,z = 0), which is taken as a representative 
value in the near field. 

and vacuum-PEC interfaces are shown, under the as- 
sumption that they are launched by slits with the same 
amplitude of the electric field on their exit apertures. 
The spatial dependence of the field is quite different in 
these two cases. Instead of the cylindrical wave- type 
algebraic decay ex \l\fx along the PEC surface, the 
field along the vacuum-gold interface shows two dif- 
ferent behaviors, separated at the crossover distance, 
Xc. It must be stressed that the chosen wavelength, 
A = 540 nm, and the type of the metal (gold) does not 
represent a special case with some particular proper- 
ties, but corresponds to a small value of Xc. Other- 
wise, the physics for this case is as rich as for other 
wavelengths. 

As shown in Fig. [3j in a very close vicinity to the 
slit X <C 1, the behavior of the field is complex and 
contains the contribution from all the angular spec- 
trum of the Green's function, or in other words, from 



all the density of EM states. Phenomenologically, the 
field in this region has been represented by a SPP plus 
an additional contribution (defined as the total field 
minus the SPP one) denoted either "creeping wave" 
or "quasi-cylindrical wave" (CW, see Refs. I5|17|) . As 
mentioned before, Eq. ([5| faithfully represents the field 
(and thus the CW) for X > 0.1 A. As X increases, 
all the smooth parts in the angular spectra in the in- 
tegral are progressively canceled out due to integra- 
tion with the oscillatory factor ~ ex.p{iqxx)^ so that 
only the sharp regions (with width Aq ~ 1/x) of the 
spectrum give a finite net contribution. These regions 
correspond to the vicinities of either the pole q = Qp 
(which is a feature of a finite width) or the kink q = 1 
(which has zero width in g-space). The field in the re- 
gion where these two contributions dominate, can be 
found by asymptotically expanding Eq. ([5|, see Sub- 
section III CI 

The electric field corresponding to the SPP and NW 
terms are rendered in Fig.jS] At distances from the slit 
of order of one wavelength, both field amplitude and 
phase [see inset (a) in Fig. [s] are well approximated 
by the SPP contribution, and the influence of the NW 
is weak. The field is locally enhanced comparing to 
the PEC case, and the efficiency of the SPP excitation 
depends upon wavelength. In the region close to the 
asymptote in the SPP dispersion relation (e^ — 1), 
the density of electromagnetic states increases and so 
it does the local field enhancement. However, while the 
mode becomes both slower and more confined, due to 
the increase of its wavevector, the absorption increases 
as well and the SPP mode is quenched at a smaller 
distance from the source. 

At distances large enough so that, due to the ab- 
sorption, the SPP is sufficiently damped, the contri- 
butions from SPP and NW are comparable, see inset 

(b) in Fig.js] This typically occurs dit Xc ^ 6 — 9Lspp^ 
with Lspp being the SPP propagation length. In the 
vicinity of Xc, the SPP and NW fields have similar 
amplitudes, so the modulus of their sum presents an 
interference behavior, leading to a set of maxima and 
minima. Notice that, in the optical region the SPP 
wavevector is close to the light cone, thus close to the 
NW one. However, qp largely increases close frequen- 
cies such that Sjri — (which, for good metals, occur 
at ujp/\/2^ where uOp is the plasma frequency), in which 
case the total field given by Eq. ([9| presents a fanciful 
two-scaled oscillatory behavior. 

The third region is located beyond Xc, where the 
contribution from the SPP field is negligible (see inset 

(c) in Fig. [3|. The field then reaches its asymptotic 
behavior, which is given by the NW, oscillating with 
the spatial period given by the free-space wavelength 
and decaying algebraically as ~ Xjx^l^ . 

Let us now analyze the dependency of the fields upon 




FIG. 4: (color online). The dependency on distance from 
the slit of the electric field modulus, at the vacuum- gold 
surface for different values of qy . The wavelength is A = 540 
nm. The normalization of the fields is the same as in Fig. [3] 
The inset (a) shows the dependency of the ratio between 
the crossover distance and SPP propagation length ratio, 
Xc/Lspp^ as a function of qy and for different wavelengths. 
The continuous curves are for Xc(qy)/Lspp{qy — 0), while 
the discontinuous ones are for Xc{qy)/Lspp{qy). The in- 
set (b) renders the normalized SPP propagation length, 
Lspp{qy) / Lspp{qy — 0) and the NW field modulus Enw^ 
together with its p- and s-polarization components, all of 
them normalized to ENw{qy — 0). 



the ^/-component of the wavevector. Fig. Q shows 
the spatial dependencies of the fields in the direction 
perpendicular to the slit for different qy. There are 
two tendencies with the increase of qy'. the amplitude 
of the field decreases, and the crossover distance Xc 
diminishes. To explain this behavior, we represent Xc 
as a function of qy (inset (a) to Fig. [4|. As can be 
seen from the curves where Xc is normalized to the 
constant ((7— independent) value Lspp{qy = 0), the 
distance after which the NW dominates decreases as 
qy increases. However, the value of Xc{qy) /Lspp{qy) 
increases as qy increases, meaning that the crossover 
occurs at smaller absolute distances, but at larger SPP 
propagation lengths, specially close to = 1. 

The inset (b) of Fig. [s] shows the NW amplitude as 
a function of qy. This behavior is due to the depen- 
dence with qy of the space-independent prefactors in 
the amplitude of the NW (which for \em\ ^ 1, goes as 
-^NW ^ ^0^^' thus decreases as qy increases). By 
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contrast, in the case of a SPP this prefactor is practi- 
cahy independent of Qy. However, the SPP propaga- 
tion length scales as Lspp oc go, so the SPP decays 
faster for larger values of Qy. As a result, as Qy in- 
creases, the NW overtakes the SPP closer to the slit, 
but with a smaller amplitude. 

As follows from Eqs.([7|),(pT]), the s-polarization com- 
ponent of the NW has a non-monotonic dependency 

upon Qy, E^^y^ oc Qy^ = Qy^e - q^. Nonetheless, 

the contribution of this component is always much 
smaller than that from p-polarization (see inset (b) of 
Fig.[4|. Notice that, as the NW presents the same alge- 
braic decay for all Qy {E^w the normaliza- 
tion to E^wiQy = 0) makes the quantities represented 
in the inset (b) of Fig. |4] independent upon distance. 
Thus, even for oblique incidence, NWs are virtually 
p-polarized waves along the interface, i.e., present the 
same polarization as SPPs. 




X/X 

FIG. 5: (color online). The dependency of the ratio Ex/Ey 
upon distance from the slit for different values of qy. The 
wavelength is A = 540 nm and the considered metal is gold. 

It is also interesting to study the ratio E^/Ey^ i.e. 
the polarization of the tangential-to-the-interface com- 
ponent of the field. For the case of in-plane launching 
{qy = 0), £^2/ = 0, and the field in the plane of the 
metal surface has only x-component. When the excited 
waveguide mode gains non-zero momentum in the di- 
rection along slit, the scattered electric field possesses a 
finite ^/-component in both far- and near-field regions. 
Fig. Is] illustrates the spatial dependency of \Ex/Ey\ 
for different qy^ showing that this ratio increases with 



qy. According to Eqs. ([6|, for a SPP the ratio E^/Ey 
scales as y^{qp/qy)'^ — 1. For a NW, from Eq. ^ 
and neglecting the contribution from the s-polarization 
(which scales as 0(1/ Sm))^ it follows that this ratio 

scales in a similar way: Ex/Ey \J'^/Qy ~ 1 • Notice 
that crossover distance is slightly different for x and 
^/-components of the electric field and, also, that the 
amplitudes of E^ and Ey are different at the crossover. 
For these reasons, the curves in Fig. |5] present dips at 
distances close to the corresponding crossovers, where 
both Ex and Ey are strongly suppressed, but Ey dom- 
inates. 



IV. CONCLUSIONS 

To conclude, in this paper we have presented the 
asymptotic expression for the EM field along the metal 
interface, launched by a subwavelength slit in the con- 
ical mount. These expressions are very accurate for 
the region down to a tenth of a wavelength. The field 
presents different contributions, which can be assigned 
to specific regions in the angular spectrum. The s- 
polarization component, although present for oblique 
incidence, not infiuence significantly the behavior of 
the field. We have studied the dependencies of the 
fields upon angle of incidence (i.e., component of the 
EM wavevector along the slit axis). In the far field re- 
gion, the distance from the slit at which the algebraic 
behavior of the field overcomes the exponential decay 
decreases as qy increases. This could favor experimen- 
tal studies of the Norton waves. 
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Appendix: Mathematical treatment of the field 
and Green's Dyadic 

In this Appendix we present the details of the ana- 
lytical computations for obtaining the field along the 
metal interface. The coordinate system is shown in 
Fig. [T] with Z = corresponding to the exit interface. 

The general self-consistent form of the field is given 
by the Lippmann-Schwinger integral equation^^ 

E(R) =Eo(R) + A:^ / dR'A£(RO^(R, R0E(R0, 
Jv 

(A.l) 
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where Eq is the solution without the sht, A£:(R) = 
£ — £m is the variation of the dielectric permittivity 
in the volume occupied by the slit V (As = outside 
the slit). The film is supposed to be optically thick, 
so that Eq = in the region of transmission, and the 
Green's function dyadic Q can be approximated by the 
one for the metal- vacuum interface, which satisfies the 
equation 

V X V X ^(R, RO - klegiYi, R^ = /(^(R - R'), 

(A.2) 

with standard boundary conditions at z = 0. In 
Eq. (A.2) / is a diagonal unit matrix; e = £ inside 
the dielectric and e = £^ inside the metal. For con- 
venience, let us introduce the dimensionless coordi- 
nates x^y^z = k^X^k^Y^k^Z. Then we assume that 
the slit is thin enough so that the field inside it has 
predominantly x-component and therefore only the x- 
column of Q will be essential, Q = e^G. This vec- 
tor is represented by a two-dimensional integral in k- 
space (see Ref. 3), and can be decomposed in p- and 
5-polarization contributions G = -\- ^ where 



ikijj 
8^ 



(A.3) 



q (Qx^Qy) with Qx^y = kx^y/kuj, Qz = \/e - , q^m 
\J £ra — and p — {x — x' — y'). The vectors a*, 



are defined to be 



qIqz 

QxQyQz 



q^qzm 



ql 

-qxqy 



(A.4) 

with t^, being the Fresnel coefficients given by 
Eq. (|3|. Taking into account that the field inside 
the slit is given by the fundamental mode and tak- 
ing into account the momentum conservation along y^ 
we can write ^{x' ^y' ^ z') = exE{x' ^ z')e^^y^y , where 
qyQ is the dimensionless ^-component of the incident 
wave wavevector. Extracting the ^/-dependency of the 
Green's dyadic, we have the following integral 



fly'^i{%o-ciy)y'+iqyy ^ 27rS{qyo - qy)e'^''y. (A.5) 



Then the integration in Eq. ( |A.3[ ) in qy is performed 
trivially. From Eq. (A.3) it follows that the integrand 
contains the exponential factor e~^^^'^^ , which decays 
at the distance of a skin depth 5 = l/lm{qzm)i and is 
of the order of a few tens of nm in the optical regime. 
We can, therefore, extend the integration limits in z' 



to [— oo, 0]. Additionally, the variation of the dyadic is 
much faster than that of the field inside the slit, hence 
the electric field inside the slit can be approximated by 
its value at the distance z = —5 (the average distance 
to the surface, weighted by the exponential decay of 
the field) 



qzm 



— £ 

(A.6) 



Then Eq. (A.l) becomes 
E{x,y,z) = C J dxG{ 



where C = 



x\y,z)Ex{x\z' = -S), 
(A.7) 

£/k^ and G = + with 
iku 



G'^P{x-x\y,z) = ^ / d^^a^'^e^[^-(^-^' 
47r J 



)+(iyy+qzz\ 
(A.8 



For brevity we have omitted "0" in qy^. We have thus 
recovered Eqs. ([T])-([3|. 

Taking into account the presence of the poles (placed 
at qz^m + qzm^ = ) and branch cuts and branch 
points (defined by l m{qz) =0), an asymptotic analysis 
of the integral (A.8) can be made following the general 



recipes described Ref. [191 as was done for the case qy 
in Ref. , Concentrating on the field at the interface 
z = 0, we simplify G in the following way 



G{x,z = 0) 



ik^ 
2^ 



r{x). 



(A.9) 



The branch cuts lm{qz) = can be removed by chang- 
ing to polar variables: qz = qo cos 0, or q^ = q^ sin (j) 
with = £ — q'^. Here we assume that the contribution 
of the branch cuts lm{qzm) = is negligible, being of 
order ~ e"'^^'^ (otherwise, some modifications in the 
solution scheme would be necessary). Then, to provide 
the exponential decay of the integrand, the variable ^ 
is further transformed into the variable s as follows: 
sincj) = 1 + zs^, so that the saddle point is placed at 
5 = 0. With this change the vector g reads 



g(x) = 
*(.) 



ds^{s)e 



d4 ^ ^[qx{s)] 
ds ' 2 



(A.IO) 



qo cos[(/)(5)]. 



where the integration path L corresponds to the real 
axis in the complex plane q. Then the singular part of 
the integrand is separated 



*o(s), *o(s) 



*(s)(s 



(A.ll) 
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where Sp = e ^g^p/^o - 1; Qxp being the x- 
component of the SPP wavevector, which satisfies, 
ql qlp qI = ££m/{£ £m)- We would hke to 
remind that only p-polarization type of poles are phys- 
ical (i.e. satisfy the radiation condition lm{qz) > 0), 
and therefore the separation (A.ll) has sense only for 



p-polarization component of The elements of the 
vector Cp are given by the residues of the integrand 
defined by Eq. 

Transforming the integration contour L to the real 
axis in the complex plane 5, the singular part of g can 
be represented using the complementary error function 
erfc 

g{x) = in Cpe^°^^*"^p^ eifc {-is p^/%x) + go(^), 



ds ( 



*o(s). 



(A.12) 



The integral appearing in the nonsingular term of g 
can be represented in the form of an infinite sum re- 
sulting from the integration of the Tailor expansion for 
^o{s)' This series reads 



go(^) 



2^ ^^d^^l'^°7~Ta±^' (^-13) 



nGeven 



(qox)- 



where F is Gamma function. For a large range of x 
only the two first terms of this expansion are impor- 
tant. With these two terms in Eq. (A. 13 ), i.e., with 
the precision up to 0(x^/^), Eq. (A.12) transforms to 
Eqs. (§-([71). 

In the far-field region, which formally occurs for 
\sp\y/qoX 1 (although comparisons with the numer- 
ical computation of the field shows that, in practice, 
this "far-field") the asymptotic expansion of the com- 
plementary error function can be used: 



erfc(- 
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(2n)! 



Spy/TTQoX 



n=0 



_^)2n+l n!(4s2(7ox)^' 
(A.14) 



where we have taken into account that Im(sp) < 0. 
Formally, this asymptotic expansion is valid for dis- 
tances such that \sp\^yqox ^ 1. However, compar- 
isons with the exact results for the fields (obtained 
form the numerical computation of relevant integrals) 
shows that usually this condition is too restrictive and 
the asymptotic expansion is valid even for shorter dis- 
tances. 
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